Vapor pressure isotope fractionation effects in 
planetary atmospheres: application to deuterium. 
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Abstract 



The impact of the vapor pressure difference between deuterated and nondeuterated 
condensing molecules in planetary atmospheres is quantitatively assessed. This difference 
results in a loss of deuterium in the vapor phase above the condensation level. In Titan, 
Uranus and Neptune, the effect on CH3D is too subtle to alter current D/H ratio determi- 
nations. In Mars, the effect can induce a large depletion of HDO, starting about one scale 
height above the condensation level. Although the current infrared measurements of the 
D/H ratio appear to be almost unaffected, the intensity of disk-averaged millimetric HDO 
lines can be modified by about 10%. The effect is much stronger in limb sounding, and can 
be easily detected from orbiter observations. 

Key words: Atmospheres, Composition; Mars, Atmosphere; Isotopes 
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1 Introduction 



Since the first detection of CH3D in the atmosphere of Jupiter (Beer and Taylor 1973), 
the D/H ratio, and to a lesser extent other isotopic ratios, have played a major role in the 
study of the origin and evolution of planetary atmospheres. Indeed, the large mass difference 
between H and its isotope D induces large differences in both the thermodynamics and 
kinetics of various processes, such as thermal escape, chemical reactions and condensation, 
resulting in an isotopic fractionation. In turn, the observed D/H ratios act as tracers of the 
chain of events which have formed the currently observed planetary atmospheres from the 
initial interstellar ices and gases. For this reason many observational efforts have focussed 
on the precise measurement of the D/H ratios in Solar System objects (see Owen (1992) 
for a review). 

However, in many circumstances, the D/H ratio observed in a particular molecular 
species does not represent the bulk composition of the planet, if any fractionation occurs 
at the sounded atmospheric levels. This is the case in the Giant Planets, where the D/H 
observed in methane (i.e. from CH3D) generally differs from that observed in hydrogen (i.e. 
HD), because isotopic exchange occurs in chemical reactions between CH4 and H2 (Beer and 
Taylor 1973, Fegley and Prinn 1988, Lecluse et al. 1996). This effect has been extensively 
studied and estimates of the corresponding fractionation factor / are available, though still 
uncertain. 

Another possible, but much less studied effect, is the isotopic fractionation caused by 
condensation, which can take place if the D/H ratio is observed in a condensible species, 
such as H2O in Mars, or CH4 in Titan, Uranus and Neptune. This effect arises from 
some slight differences between the vapor pressures of isotopes of a given chemical species, 
with the heavy isotope molecule generally having a lower vapor pressure. During cloud 
formation in a planetary atmosphere, the condensed phase is then enriched, and the vapor 
phase depleted, in the heavy isotope. The depletion increases with the amount of condensed 
mass. This effect is widely documented on Earth, both theoretically and observationally, 
with aircraft measurements determining a D/H ratio as low as 0.7 times the standard 
ratio (SMOW) a few kilometers above the condensation level (Taylor 1972). This raises 
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concern that the D /H ratios measured in other planetary atmospheres could also suffer from 
significant departures from the global planetary ratios, due to condensation fractionation. 

The subject of this paper is to investigate the effect of the isotopic fractionation occuring 
in cloud condensation, and how it affects the measured D/H ratio. We focus on water in the 
atmosphere of Mars, and on methane in Titan, Uranus and Neptune. In the next section we 
present the theoretical bases of isotopic variations induced by cloud formation, along with 
the available laboratory data relevant to the studied planets. Then we present the expected 
isotopic variations and the consequences on the observed D/H ratios. In the case of Mars, 
we also present synthetic calculations of millimetric HDO lines in order to determine the 
detectability of this isotopic fractionation. 

2 Laboratory measurements and isotopic cloud models 

2.1 Fractionation coefficient 

The difference in equilibrium vapor pressures of isotopic molecules, often referred to as 
the vapor pressure isotope effect (VPIE), exhibits a large variety of behaviors, revealing the 
nature of the intramolecular forces at work in the condensed phases. For this reason, many 
laboratory and theoretical studies have been aimed at precisely measuring or modeling the 
VPIE, in particular for water and methane. A comparison of the different values can be 
found in Jancso and Van Hook (1974), but here, following Jouzel (1986) and Bigeleisen et 
al. (1967), we adopt the fractionation coefficients of Merlivat and Nief (1967) for water and 
of Armstrong et al. (1953) for methane. 

The isotopic content of an atmospheric parcel can be expressed in terms of 5: 



where (D/H) D is the standard atmospheric value, i.e. the global D/H ratio of the planet. 
The VPIE is expressed in terms of a fractionation coefficient a, defined as the ratio of the 
D/H in the condensed phase to that in the vapor phase: 
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where the c and v subscripts refer respectively to the condensed and vapor phases. 

The water fractionation coefficient for the solid/vapor transition is given by Merlivat 
and Nief (1967): 

16.288 o o 
lna= T2 x 10 3 - 9.34 x 10~ 2 , (3) 

where T is the temperature in Kelvin. Unfortunately, this expression is experimentally 
validated only between 230 and 360 K. No data are available at the temperatures relevant 
for Mars' atmosphere, where water condensation occurs typically between 190 K and 140 K. 
The extrapolation of the above formula down to these temperatures might be hazardous. 
Therefore, we explore the consequences of isotopic fractionation in two cases. In the first 
one, we extrapolate the expression of Merlivat and Nief (1967) to the full Martian tem- 
perature range (a = 1.43 at 190 K, and a = 2.09 at 140 K). In a second, conservative 
approach, we adopt the coefficient at 230 K (a = 1.24) for all lower temperatures. This 
likely underestimates the effects of isotopic fractionation. 

The methane fractionation coefficient has been measured for the solid/vapor transition 
by Armstrong et al. (1953): 

110.2 1.26 

lna = — 2 — , (4) 

where T is the temperature in Kelvin. This expression is based on laboratory measurements 
carried out between 75 and 91 K. The temperature minimum in Uranus and Neptune lies 
around 55 K. However at this temperature the extrapolated fractionation coefficient is still 
very low (a = 1.014). In addition, Bigeleisen et al. (1967) gave a theoretical justification to 
the expression of Armstrong et al. (1953). Therefore, we think it can safely be extrapolated 
down to 55 K. Methane exhibits the rare behavior of an inverse isotope effect: the condensed 
phase is enriched in the lightest isotope (i.e. CH4) at temperatures warmer than 87.5 K, 
while at lower temperatures the "normal" effect (i.e. enrichment in CH 3 D), takes place. 
However, as in Titan's, Uranus' and Neptune's atmospheres methane condenses only at 
lower temperatures, the global effect is a depletion of deuterium in the gas phase, as we 
show below. 
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2.2 Cloud models 

Cloud formation occurs when a parcel of moist air cools below its dewpoint. Cooling 
can be due to a variety of processes: radiative heat loss, mixing of different air masses, 
expansion due to vertical lifting. In the planetary atmospheres considered here, most of 
cloud condensation is caused by vertical lifting, either in convection for Uranus and Neptune 
(Carlson et al. 1988), or in global circulation Hadley cells on Mars (Jakosky and Haberle 
1992) and Titan (Hourdin et al. 1995). For this type of cooling, the isotopic fractionation 
theory has been extensively developed for the Earth's atmosphere (see Jouzel (1986) for a 
review). Modeling of actual clouds is difficult, due to the fact that part of the condensed 
phase is removed by precipitation, while the remaining part ascends with the updraft. 
Moreover, in planetary atmospheres the condensation processes are still poorly understood, 
compared to the Earth's atmosphere. Therefore we restrict our study to two simple ideal 
cases considered by Dansgaard (1964) and designed to span the whole range of possible 
cloud models, and thus of isotopic fractionation values. 

In the following, we assume that the condensed and vapor phases are always in isotopic 
equilibrium, i.e. that the isotopic equilibrium is kinetically achieved. This is clearly not 
the case in the Earth's atmosphere (Taylor 1972). Jouzel et al. (1980) have shown that 
neglecting kinetics effects, in the case of hailstone formation in the Earth's atmosphere, 
leads to a 20% overestimation of the D/H depletion. However, modeling this departure is 
beyond the scope of this work which aims at a first order estimate of the VP IE in planetary 
atmospheres. 

The first model is based on an open cloud system (Dansgaard 1964), where the solid 
phase is assumed to condense in isotopic equilibrium with the vapor phase, and to leave the 
cloud by precipitation immediately after its formation. In this case, the elementary isotopic 
variation of the vapor phase is given by: 





leading at first order to: 
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where n v and n l v are the vapor mixing ratios respectively of the nondeuterated and deuter- 
ated species, and dn v and dn l v their elementary variations. Since the condensed phase is in 
isotopic equilibrium with the vapor phase we have the relation: 

^k = a Vk (7) 

dn v n v ' 



from which it results that: 



\ a ~ 1 ) — - ( 8 ) 



1 + 5 V ' n v 

The open cloud system gives the largest possible depletion of the gas phase in deuterium 
above the condensation level, since the isotopic excess in the condensed phase is definitely 
and constantly removed by precipitation. 

In contrast, the second model induces the smallest possible depletion. It is based on 
a closed cloud system (Dansgaard 1964), where the condensed and vapor phases ascend 
together in the air parcel. In this case the number of molecules is conserved: 

( n Q = n c + n v for the nondeuterated species 

\ n 5 = n v 5 v + n c 5 c for the deuterated species 

where n v (respectively n c ) is the gas mixing ratio of the nondeuterated species in the 
vapor phase (respectively condensed) phase, and n and 5 are the mixing ratio of the 
nondeuterated species and the isotope ratio below the condensation level. It thus results 
from the combination of Eq. || (isotopic equilibrium) and Eq. [J that the isotopic ratio of 
the vapor phase is given by: 

n (5 - a + 1) + (a - l)n v , . 

°v = 7Z \ ; • (10) 

(1 — a)n v + an a 

These two theoretical models have been experimentally validated for the terrestrial 
atmosphere. Taylor (1972) has shown from in-situ measurements that the deuterium content 
of air masses above the condensation level was in agreement with the open cloud model, 
while laboratory analyses of hailstones (Jouzel et al, 1975) support the closed cloud model 
under stormy conditions. In addition, we have validated our code by running it in terrestrial 
conditions, and comparing our results with the calculations of Jouzel (1986) (their Fig. 4). 
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3 Atmospheric modeling 
3.1 Mars: HDO 

The Martian D/H ratio was determined in the infrared by Owen et al. (1988) and 
Krasnopolsky et al. (1997) at 3.7 /um, and by Bjoraker et al. (1989) at 2.65 /im. In this 
solar reflected part of the spectrum, the observations are sensitive to the whole column 
density, and the D/H ratio is estimated by comparing the total column densities of H2O 
and HDO. For this reason, the impact of the VPIE on the D/H determination strongly 
varies with the altitude of water condensation. If condensation occurs close to the surface a 
large fraction of the HDO column can be affected by the isotopic fractionation. If it occurs 
a few scale heights above the surface the modification to the column is only marginal. A 
proper estimate of the VPIE on Mars must then take into account the variations of the 
water condensation level, due to variations in temperatures and water content of the lower 
atmosphere. 

The Martian climate exhibits strong annual variations due to the high eccentricity of 
the planet's orbit, on which are superimposed episodes of global warming correlated with 
the dust content of the atmosphere (Clancy et al. 1996). We investigate the effect of the 
isotopic fractionation in the case of three different temperature profiles. In the lower at- 
mosphere, we select profiles from the disk-averaged measurements of Clancy et al. (1996). 
Such observations are most appropriate to model the effects of isotopic fractionation as 
they would be observed from ground based observations of the planet. The millimeter 
wave observations of Clancy et al. (1996) typically constrain the atmospheric tempera- 
ture up to 45 km (0.05 mbar). At higher levels, we take temperature information from 
the numerical simulations of the Martian Climate Database (MCD) (Read et al. 1997, 
http://www.lmd.jussieu.fr/mars.htmi ). Our three thermal profiles are the following: 



Cold profile: at Mars' aphelion, in the case of a low dust loading of the 
atmosphere, we adopt the coldest temperatures of Clancy et al. (1996): 
190 K at 3 mbar, 160 K at 0.5 mbar and 140 K at 0.05 mbar (Fig. la). 
Following the MCD, the temperature is left constant at 140 K above this 
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pressure level. 

• Mean profile: at Mars' perihelion, in the case of a low dust loading of the 
atmosphere, the temperatures are 210 K at 3mbar and 180 K at 0.5mbar 
(Fig. 2a). Following the MCD, we place the tropopause at 0.05 mbar with 
a temperature of 140 K. 

• Hot profile: at Mars' perihelion and with a high dust loading, we adopt 
a temperature profile in agreement with the April 1992 measurement of 
Clancy et al. (1996) and with that observed by Viking/IRTM (Martin 
and Kieffer 1979). This profile has 220 K at 3 mbar, 200 K at 0.5 mbar, 
and 180 K at 0.05 mbar. The tropopause is fixed at 5/ubar and 140 K. 

The amount of water in the Martian atmosphere is highly variable, both spatially and 
seasonally (Jakosky and Haberle 1992), as a result of exchanges between the atmosphere 
and non-atmospheric reservoirs, such as polar caps or ground ice. In order to investigate the 
largest set of possible situations we calculate the isotopic fractionation for each of the three 
temperature profiles with three typical water column densities: 3, 10 and 30 precipitable 
micrometers (pr./im). 

Although atmospheric water vapor is in equilibrium with the ice reservoirs, the global 
(i.e. below the clouds) atmospheric D/H ratio is not expected to differ from that in water ice. 
This is because fractionation effects occur at each phase change except at sublimation and 
melting of ice (Jouzel 1986). This behavior comes from the very slow molecular diffusion 
within ice, which makes the time scale of isotope exchanges inside the ice longer than the 
time scale of phase change (Joussaume et al, 1984). (For the same reason, the D/H ratio 
measured in cometary comae should represent the true nuclear value.) 

For the three different temperature profiles and water column densities we calculate the 
expected vertical profiles of the D/H ratio (more precisely of S v ). The results are displayed 
on Fig. 1, Fig. 2 and Fig. 3. In the modeling, we assume that the condensation process 
occurs at 100% relative humidity, i.e. we do not allow for supersaturation. The closed cloud 
model gives rise to a moderate depletion in HDO. <5^ asymptotically converges toward —19% 
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at the top of the atmosphere if a is kept constant to its value at 230 K. In the temperature- 
dependent case the depletion is larger: 5 V = —52%. These values can be inferred from 
the asymptotic value of Eq. |l(]: 6 v = (1 — a) /a, where a is taken at 230 K and 140 K 
respectively. In contrast, the open cloud model induces a huge depletion in HDO in both 
cases. The D/H ratio at the tropopause is ~ 20% of its value at the surface (5 V ~ —80%) 
if a is kept constant, and only ~ 1% (5 V ~ —99%) in the temperature-dependent case. 

To estimate the effect of the isotopic fractionation on the D/H determinations, we 
compare the HDO column density from the models with that expected for an altitude 
independent D/H ratio. The results of this calculation are given in Table I. As expected, 
the column integrated HDO depletion strongly depends on the level of cloud formation, and 
consequently on the temperature profile. The hot and mean profiles lead to a very small 
depletion, less than 1% (Table I). On the contrary, a depletion of the HDO column of the 
order of 10% occurs for the cold profile, i.e. when the condensation level lies close to the 
surface. Note that the depletion in the HDO column is less sensitive to the H2O amount 
than to the temperature profile. 

Because of the uncertainty on the behavior of real martian clouds and on the extrap- 
olation of the fractionation coefficient, it seems difficult to accurately correct for the con- 
densation effects. However, in what follows, we examine to which extent past and future 
observations of HDO/H2O may be affected. 

Measurements of the Martian D/H ratio have been reported by Owen et al. (1988) and 
Krasnopolsky et al. (1997). However, their observational uncertainties, D/H = 6 ± 3 and 
5.5 ± 2 times the standard terrestrial ratio (SMOW), are so large that taking into account 
the fractionation effect has little value. Bjoraker et al. (1989) obtained observations of Mars 
near perihelion (August 5, 1988), and reported a much more precise value of 5.2 ±0.2 times 
the standard terrestrial ratio. Since these data have not been published extensively, it is 
difficult to assess their accuracy. In any case, at the time of these observations, the observed 
water abundance was close to 10 pr. /im, and the temperature measured at 5.1mbar was 
high (228 K). Our hot profile is thus the most appropriate to estimate the HDO column 
depletion, showing that the reported D /H ratio is essentially unaltered by the fractionation 
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effect. 

HDO has been also detected in the microwave range (Encrenaz et al. 1991, Clancy 
et al. 1992), and used as an estimate of H2O assuming the D/H ratio of Bjoraker et al. 
(1989). The vapor pressure isotope effect could bias the water abundance determination. 
Comparison of HDO and H2O or H2 18 observations in this spectral range could also offer 
an opportunity to detect this effect. This is shown in Fig. 4, where we compare synthetic 
profiles of the 226 GHz HDO line, calculated with an altitude independent D/H profile, 
and with the D/H profiles generated in the case of the cold thermal profile and a water 
abundance of 10 pr. fim. In this calculation, we assume a whole-disk observation. For this 
optically thin line (r ~ 10 -2 ) the emission mostly originates from the limb of the planet. 
We find that the HDO depletion induces a significant change in the line intensity. The open 
cloud model with a temperature-dependent fractionation coefficient decreases the contrast 
of the HDO-line by ~ 15%. The other three models decrease the contrast by 5-10%. Such 
an effect is at the limit of the detectability of current observations. 

A much more important effect arises if one exclusively probes atmospheric levels above 
the condensation level from orbiter observations in horizontal viewing. We simulate such 
limb observations in several cases spanning the different temperature profiles and water 
abundances. For this, we compare the HDO emission of two D/H profiles: a vertically 
independent D/H ratio, and that predicted by the closed cloud model with a temperature- 
dependent a (since it represents the average of our four cloud systems). The synthetic 
calculations are displayed on Fig. 5. They clearly show that the modeled D/H depletion 
leads to important differences in the line intensity from the constant D/H ratio case. The 
reason is that, one scale height above the condensation level, already 20 to 30% of the 
deuterium is lost. Thus, probing at this altitude induces a 20 to 30% decrease in the HDO 
line intensity. The effect is even larger if the tangent altitude is higher. 

There are current plans to fly a millimeter heterodyne instrument in Martian orbit, 
particularly in the framework of Martian micromissions. Our calculations show that if 
HDO lines are to be used to retrieve the vertical profile of H2O, the vapor pressure isotope 
effect needs to be considered. Conversely, comparing limb measurements on HDO and H2O 
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(or H2 18 0) is a promising strategy to detect and measure this effect. 

3.2 Neptune, Uranus and Titan: CH 3 D 

The D/H ratio in Neptune has been measured from CH3D both in the troposphere at 
1.6 /im (de Bergh et al. 1990) and in the stratosphere at 8.6 ^m (Orton et al. 1992, Bezard 
et al. 1997). The results of these three studies are compatible within their error bars, de 
Bergh et al. (1990) obtain (9^ 2 ) x 10~ 5 while Orton et al. (1992) give (5.6±J;|) x 10~ 5 and 
Bezard et al. (1997) (5t£f) x 10~ 5 . (Here, the values of de Bergh et al. (1990) and Orton 
et al. (1992) have been updated to account for the new value of the isotopic enrichment 
factor between H2 and CH4 recommended by Lecluse et al. (1996), / = 1.61 ± 0.21, cf. 
Feuchtgruber et al. (1999); we note that this value of / is significantly different from that 
proposed by Fegley and Prinn (1988)). These three determinations can be affected by the 
VPIE, since the CH^D-lines probe atmospheric levels well above the tropopause at 8.6 /jm 
or around the condensation level at 1.6 ^m, although the effect must be much smaller in 
the latter case. In Uranus the only detection of CH3D was obtained at 1.6 ^m (de Bergh 
et al. 1986), yielding D/H = (5.4±|;g) x 10~ 5 for / = 1.68 ± 0.23 (Lecluse et al. 1996, 
Feuchtgruber et al. 1999). Similarly to Neptune, the sounded atmospheric levels lie near 
the methane condensation level, and the D/H ratio can, in principle, be affected by isotopic 
fractionation. 

Testing for this possibility, we model the CH4 condensation process in Neptune using the 
temperature profile of Feuchtgruber et al. (1999) (Fig. 6a). We adopt the deep methane 
mole fraction of Baines et al. (1995), /ch 4 = 0.022. As for Mars, we assume that ice 
formation occurs at 100% relative humidity. Figures 6b and 6c respectively display the 
fractionation coefficient a and the calculated methane vertical distribution as a function of 
pressure. It can be seen that methane condensation occurs at too shallow a level (1.5 bar, 
80 K) to give rise to an inverse isotope effect. Thus, as for Mars, a loss of deuterium in the 
gas phase results is the net result, as shown on Fig. 6d. 

According to Carlson et al. (1988) the rapid condensation of methane implies massive 
precipitation. In this case the open cloud system is the best model of the actual D/H 
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depletion. Yet, the stratospheric measurement of methane, /ch 4 = (l-^lo!?) x 10~ 3 (Bezard 
1998), indicates a supersaturation with respect to the tropopause cold trap value (1.3 x 
10 -4 ), implying that some ice particles are transported up to the stratosphere where they 
sublimate. We therefore adopt as a stratospheric estimate of the D/H depletion the value 
obtained at the pressure level, 0.7 bar, where the methane mole fraction is equal to 1.3x 10 -3 . 
This yields 5 V = —1% (Fig. 6d). Even if the observed supersaturation of CH4 is not 
accounted for, the asymptotic (stratospheric) value of 5 V is only —2.8% (Fig. 6d). Thus 
the effect is too subtle to alter the interpretation of current measurements at 8.6 fim. This 
is naturally even more true for the 1.6 /im measurements. The same conclusion holds for 
Uranus, since its temperature profile and its methane abundance are similar to those of 
Neptune. 

In Uranus and Neptune, the D/H ratio is also measured from HD (see Feuchtgruber et 
al. (1999)). However, comparing the measurements from HD with those from CH3D can- 
not be used presently as an observational determination of the condensation fractionation 
effect, not only because the measurements are not accurate enough, but also because the 
fractionation factor between CH4 and H2 is not known with enough accuracy (Fegley and 
Prinn 1988, Lecluse et al. 1996). 

In Titan there is no definite evidence for a permanent methane cloud, although Griffith 
et al. (1998) have reported observations of transient clouds near 15 km. The methane 
tropospheric abundance is also not yet accurately determined. Many different methane 
profiles have been derived (Flasar 1983, Lellouch et al. 1989, Courtin et al. 1995). McKay 
et al. (1997) were even able to fit the Voyager occultation data with a constant CH4 mixing 
ratio throughout the whole atmosphere. In this situation, we only calculate an upper limit 
of the VPIE in Titan. 

For this, we assume that the fractionation coefficient is constant in the whole troposphere 
at its highest possible value, i.e. its value at the tropopause. Adopting the temperature 
profile of Yelle et al. (1997) with a tropopause at 71 K, we find a = 1.004. Within this 
assumption Eq. || can be solved analytically: S v = (1 + $o){n v /no) a ~ l — 1. Flasar (1983) 
pointed out that the surface humidity must be lower than 0.7 to keep the troposphere 
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convectively stable, but here we adopt the conservative approach of a methane profile 
saturated throughout the troposphere. This again can only enhance the effect, and yields 
a ratio of the mixing ratio at the surface to that at the tropopause of n v /no = 0.21. All 
these assumptions result in an upper limit of S v = —0.6% to the D/H depletion in the 
stratosphere, totally negligible for the D/H determination from stratospheric methane. 

4 Conclusions 

We have investigated the consequences of the vapor pressure isotopic effect on con- 
densible deuterated species in planetary atmospheres. This effect induces a depletion of 
deuterium in the vapor phase above the condensation level. The magnitude of the deple- 
tion largely depends on the type of cloud system. In the atmospheres of Titan, Uranus and 
Neptune, the depletion in CH3D is not detectable given current observational uncertainties, 
so that D/H ratios determinations are not altered by this isotopic fractionation. 

In the atmosphere of Mars HDO is greatly affected. Depletions from 20% to 99% of the 
D/H ratio are possible a few scale heights above the condensation level. The magnitude of 
the effect depends primarily on the temperature profile. In particular the D/H measurement 
of Bjoraker et al. (1989), performed during a relatively warm period, does not suffer from 
isotopic fractionation. When the condensation level lies close to the surface, the isotopic 
fractionation effect induces variations of ~ 10% in the HDO-line intensity for disk averaged 
millimetric observations. Such an effect might become detectable in the near future as the 
accuracy of these observations will improve. 

The effect is much more spectacular in horizontal viewing. Limb sounding above the 
condensation level thus appears to be an efficient method to detect the HDO depletion. It 
is an interesting scientific objective for future Martian missions. Indeed, if detected and 
quantified, the vapor pressure isotope effect could enhance our understanding of the water 
condensation process and of the atmospheric circulation. 
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Table I: Loss of HDO in Mars' atmosphere: difference between the VPIE-fractionated 
column densities and the non-fractionated column density. 



Water Open Model Closed Model 

abundance (pr. /an) Cold Mean Hot Cold Mean Hot 

Constant 3 1.3% 0.19% < 0.1% 0.81% 0.11% < 0.1% 

a 10 3.1% 0.33% < 0.1% 1.9% 0.20% < 0.1% 

30 10% 0.65% < 0.1% 6.4% 0.40% < 0.1% 

Temperature 3 2.5% 0.38% < 0.1% 

dependent 10 5.1% 0.63% < 0.1% 

a 30 15% 1.1% < 0.1% 



1.6% 0.25% < 0.1% 
3.4% 0.41% < 0.1% 
10% 0.73% < 0.1% 
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Figure Captions 



Figure 1 D/H depletion in Mars' atmosphere for the cold temperature profile. Panel a: 
Temperature profile. Panel b: S v vertical profiles calculated for a water abundance of 3 
pr. fim with an open cloud model for temperature-dependent a (0M1, solid thick line) and 
for a constant a (0M2, dashed thick line), with a closed cloud model for temperature- 
dependent a (CM1, solid thin line) and for a constant a (CM2, dashed thin line). Panels 
c and d: Same for 10 and 30 pr. fim. 

Figure 2 Same as figure 1 for the mean temperature profile. 
Figure 3 Same as figure 1 for the hot temperature profile. 

Figure 4 Comparison of disk averaged 226-GHz HDO line in Mars' atmosphere for a non- 
fractionated profile (tripledotted-dashed line) with that obtained for the four different HDO 
profiles calculated with the cold profile and 10 pr. /xm(see Fig. lc): OM1 (solid thick line), 
OM2 (dashed thick line), CM1 (solid thin line) and CM2 (dashed thin line). 

Figure 5 Comparison of limb observations of the 226-GHz HDO line calculated with a 
non-fractionated profile (solid line) with that calculated with a closed cloud model and a 
temperature-dependent a (CM1, dashed-dotted line) in four different cases. Panel a: Cold 
temperature profile, water abundance of 30 pr. /xm and a tangent altitude of 5 km. Panel b: 
Cold temperature profile, water abundance of 10 pr. /j,m and a tangent altitude of 20 km. 
Panel c: Mean temperature profile, water abundance of 10 pr. /jm and a tangent altitude 
of 30 km. Panel d: Cold temperature profile, water abundance of 3 pr. fim and a tangent 
altitude of 60 km. 

Figure 6 D/H depletion in Neptune's atmosphere. Panel a: Temperature profile. Panel 
b: Fractionation coefficient as a function of pressure. Panel c: CH4 vertical profile. Panel 
d: D/H vertical profile calculated with an open cloud model (solid line) and a closed cloud 
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model (dashed line). 
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